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We briefly review a number of important recent experimental and theoretical developments in 
the held of dynamic fracture. Topics include experimental validation of the equations of motion for 
straight tensile cracks (in both inhnite media and strip geometries), validation of a new theoretical 
description of the near-tip helds of dynamic cracks incorporating weak elastic nonlinearities, a new 
understanding of dynamic instabilities of tensile cracks in both 2D and 3D, crack front dynamics, and 
the relation between frictional motion and dynamic shear cracks. Related future research directions 
are briehy discussed. 


I. INTRODUCTION 


The goal of this paper is to offer some perspectives on a number of signihcant developments in the held of 
dynamic fracture that have taken place over the past decade or so. The subjects that we chose to describe are 
topics that we have worked on, so this paper is, by dehnition, quite subjective. We believe that in these chosen 
areas there has been signihcant progress — and each is interesting in its own right. We will make no attempt to 
present any of these subjects in depth. The idea is to present a schematic, non-technical, overview of the main 
discoveries and fundamental insights that have evolved from recent research in each of these different areas. 
Many of these topics have been driven forward by new experiments that have utilized a variety of new techniques 
to obtain measurements that had not previously been possible. In these cases, comparison of these experiments 
to the existing theoretical predictions revealed previously unappreciated gaps in our understanding. These, in 
turn, provided incentives to expand the theoretical framework to encompass these observations. Some, but not 
all, of the topics addressed in this perspective piece were systematically reviewed in two relatively recent review 
papers mm- Our goal here, as explained above, is different. 

The structure of this paper is as follows. In Sect. [H] we consider the equation of motion for dynamic tensile 
(mode I) cracks that has been predicted in the framework of Linear Elastic Eracture Mechanics (LEEM). We 
will describe recent experiments that have quantitatively shown that the ideas of energy balance that had been 
the cornerstone of dynamic fracture mechanics indeed provide an excellent quantitative description of crack 
dynamics as long as rapid fracture is mediated by single “simple” cracks; i.e. cracks that propagate along 
prescribed paths and have not undergone dynamic instabilities. 

Section [ni] considers the form of the singular fields surrounding the tips of these “simple” cracks. We first dis¬ 
cuss a few rather peculiar features of the ^/r singular fields predicted by LEEM. We then show that experiments 
which utilize soft materials to “slow down” crack dynamics, making them accessible to direct measurements, 
revealed that the fields in the vicinity of dynamic crack tips may significantly deviate from the LEEM pre¬ 
dictions. These fields were explained by theoretically incorporating the leading nonlinear elastic behavior of 
the medium that becomes important at the high strains that invalidate LEEM as the tip is approached. This 
theory also predicts the existence of a new intrinsic length scale emerging from the competition between linear 
and nonlinear elastic constitutive behaviors. Thus, while energy balance predicts the overall dynamics of simple 
cracks, the structure of the elastic fields driving fracture is modified by elastic nonlinearities near the crack’s 
tip, where linear elasticity starts to break down. 


In Sect. 1^ we go on to consider tensile cracks that are no longer simple. This happens when simple cracks 
become unstable to different instabilities. We briefly mention the micro-branching instability, where frustrated 
microscopic branches are spawned at the crack tip, highlighting both its role in significantly limiting crack 
velocities and its 3D nature. We then show that when micro-branching is suppressed by approaching the 
2D limit, cracks accelerate to extremely high velocities, where they become unstable to spontaneous shape 
oscillations at extreme velocities. Recent theoretical and experimental developments have shed new light on 
the origin of this instability, and its relation to crack tip nonlinearities through an intrinsic length scale. 

The results discussed in Sects. [n]|IV| mainly focus on cracks that propagate in an effectively 2D medium. 
Section |V| discusses phenomena that entail an intrinsically 3D description of fracture dynamics. Based on 
very recent experiments, we provide additional insight into the properties and origin of the micro-branching 
instability and provide evidence that it is related to the oscillatory instability. We then briefly discuss both 
in-plane and out-of-plane crack front dynamics. 

While Sects. HYI consider cracks that are ostensibly tensile in nature, Sect. |VI| briefly reviews studies that 
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reveal relations between frictional motion and dynamic shear (mode II) cracks. Experiments over the last 
decade or so have revealed tight connections between interfacial friction dynamics and crack-like ruptures that 
appear similar to earthquakes along natural faults. Recent experiments, moreover, suggest that at least some 
of these can be described with classical mode II crack solutions. Understanding others requires richer coupling 
of bulk elasticity to constitutive descriptions of friction. 

We conclude with some opinions on various directions in which the field of dynamic fracture may evolve in 
the future. 


II. THE EQUATION OF MOTION FOR SIMPLE CRACKS 

What are “simple” cracks? Let us define a simple crack as a single crack that propagates along a single 
(straight) fracture plane [2]. A simple crack is, basically, a crack that does not undergo any path instabilities. 
How is the motion of such cracks described? The equation of motion for simple cracks is based on the principle 
of energy balance llE]; the energy release rate G (the elastic energy released per unit crack surface) is balanced 
by the fracture energy E (the energy dissipation per unit crack surface). 

As long as the surrounding material is elastic, F represents all of the energy dissipated by a crack. For brittle 
materials, the “small scale yielding” assumption is valid; all of the dissipative processes are encompassed in a 
sufficiently small region surrounding the crack tip [3]. T{v) is, in general, a material function that may depend 
on the crack velocity, v HHig. Its i;-dependence is often unknown, but in some cases can be measured [IMS. 

The other “leg” of the energy balance equation is the energy release rate, G. The latter is a property of 
the elastodynamic solution of a moving crack and can be obtained by the J-integral for steady-state cracks 
or more generally by the dynamic stress intensity factor. Within the framework of Linear Elastic Fracture 
Mechanics (LEFM), the near tip stress cr (as well as displacement gradient) fields are characterized by a 
universal singularity, cr ^ Ki^n /where the intensity of the singularity is quantified by the stress intensity 
factors Kijj. Kj corresponds to mode I fracture, i.e. tensile cracks, and Ku corresponds to mode II fracture, 
i.e. shear cracks. The stress intensity factors depend on the loading configuration and geometry of a given 
fracture problem. Once they are known, the energy release rate G can be readily calculated. The equation of 
motion for a straight crack is then given by T{v) = G. Next, we briefly discuss two such equations, which were 
quantitatively tested in recent experiments. 

The most well-known equation of motion was obtained for a finite crack propagating in an infinite medium. 
In this case, the equation takes the form [3] 

, ( 1 ) 

where a pre-factor of order unity was omitted. Here I is the crack length {l = v)^ (t^o is the externally imposed 
tensile stress far from the crack, E is Young’s modulus and cr is the Rayleigh wave-speed. Once r{v) is 
known, Eq. 0 is a first-order nonlinear evolution equation for l{t). It is valid as long as “small scale yielding” 
conditions prevail (i.e. linear elasticity is valid everywhere except for a small region near the tip, the so-called 
process zone). Additionally, the externally applied stresses should be able to be mapped to tractions applied to 
the crack faces and the crack cannot interact with the external boundaries (i.e. it propagates in an effectively 
unbounded medium). 

The most interesting properties of Eq. 0 are that the crack tip dynamics are not inertial, i.e. the tip behaves 
as a massless “particle” where the acceleration v plays no role (of course material inertia is essential for this 
equation to hold), and that the limiting propagation velocity is cr. The equation has been shown to beautifully 
describe crack’s motion when its underlying assumptions are met [HiiiMinj- In particular, this description 
was shown to be valid in experiments in the brittle acrylic, PMMA, for v < OAcr (when cracks are indeed 
simple) [17] . More recent experiments showed the equation of motion to be in excellent quantitative agreement 
with measurements of dynamic cracks in brittle gels for an unprecedented range of velocities, 0 <v <0.96cr, 
approaching the limiting velocity m- This wide velocity range was only accomplished by suppressing crack 
instabilities. This excellent agreement with theory, presented in Fig.[^ was obtained entirely from first principles 
with no adjustable parameters. 

Less known is an equation of motion for simple cracks which is valid when the crack “knows” about the 
external boundaries. While energy balance is always respected, Eq. 0 was derived for an system of infinite 
extent and, therefore, does not account for energy reflected back to the crack from external boundaries. When 
the waves that the crack emits during propagation are reflected back from the outer boundaries, the crack’s 
dynamics may undergo a qualitative change. Let us consider the important particular case of a crack propa¬ 
gating in a long strip of width 2b which is loaded uniaxially by fixed tensile displacements. When the crack is 
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FIG. 1: The measured crack propagation velocity i’ as a function of the crack length I (circles) compared to the LEFM 
predictions in Eq. 0 (solid lines), for various values of the far field tensile stress ctoo = 15.7,14.9,12.4,10.4kPa (top to 
bottom). All measurements are for single cracks propagating in effectively “infinite medium” conditions (/(t) <2b — 
125mm). Adapted from [T]. 


sufficiently long, I ^ 6, there is a fixed amount of strain energy per unit area W stored way ahead of the tip and 
no stored energy well behind the tip. Consequently, under steady state conditions the crack velocity is simply 
determined by F(u) = W. Marder derived the crack evolution approaching this steady state to leading order 
in the dimensionless acceleration 5u/c^, where Cd is the dilatational wave-speed [20]. In this case, the crack 
length plays no role (the relevant length scale is 6, not /), but the acceleration (which quantifies the deviation 
from steady state) is the relevant dynamical quantity. It is worth noting that this acceleration dependence is 
entirely due to the crack’s interaction with the system boundaries and not to, for example, large acceleration 
corrections to the stress fields predicted for rapidly accelerating cracks in an infinite system [21]. This (little 
known) equation of motion takes the form [20] 


r{v) = G{v,v)-W 



( 2 ) 


This equation is qualitatively different from Eq. 0 in one major aspect: the crack tip dynamics are effec¬ 
tively inertial, involving acceleration and an “effective mass” proportional to (1 rather than being 

“massless” as in Eq. 0 . These effective inertial dynamics are inherited from the interaction of the crack with 
the strip boundaries (and hence with its own history, giving rise to “memory”). Since W is externally controlled 
and F(u) is a bounded material function (attaining a finite value as v^cr), one can impose W>r{cR). In this 
case, the divergence of the effective mass proportional to (1 as V ^ cr ensures that cr is still the 

limiting velocity as in Eq. 0 (here such that v{l — I(?r) ^ remains finite as v^cr). The predictions 

of Eqs. Q-Q can be tested in a single experiment involving a long strip under fixed displacement boundary 
conditions. For short times and /(t) ^ 6, we expect Eq. 0 to be valid, while for longer times and /(t) ^ 6, 
we expect Eq. 0 to be valid, with some crossover in between. A comparison to such an experiment, shown 
in Fig. demonstrates that both Eqs. 0-0 are quantitatively correct when valid and that the transition 
between them is rather sharp. 

Therefore, while for many the equation of motion described by Eq. 0 is unjustifiably thought to be universal, 
Eq. (|§ explicitly demonstrates that this is not the case. In fact, we fully expect that the qualitative dynamics 
which are inherent in Eq. Q are in many ways general. If, for whatever reason, a crack is allowed to interact 
with its past history, we should expect strong memory, inertial-like, effects. Examples where such interactions 
may occur are when a crack tip/front encounters a localized inhomogeneity or obstacle in its path. In such 
cases, when the translation invariance along the singular crack front — the crack’s leading edge in 3D — is 
broken, it has been shown that waves will propagate along the crack’s front [22H28] . These waves, which transfer 
energy between different points along the crack front, may influence crack dynamics in much the way that the 
reflected waves that lead to Eq. 0 impart a crack with inertia. Moreover, experiments have shown that the 
crack front indeed retains a memory of an obstacle after it has passed it; damped oscillations of the front are 
observed immediately in the wake of a localized obstacle [26]. These issues will be briefly discussed in Sect. [Vj 
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FIG. 2: The measured crack velocity vs. crack length in a sample of dimensions 120 x 30 x 0.2mm (circles) compared to 
the LEFM prediction for an inhnite medium (small /(t)), Eq. 0 (red line), and to the LEEM prediction for an inhnite 
strip (large /(t)), Eq. 0 (black line). The rather sharp transition between the two regimes occurs at t^2b/cs {b is the 
half-width of the strip and Cs is the shear wave-speed), when rehected waves from the vertical boundaries interact with 
the crack tip. Adapted from HU. 


III. WEAKLY NONLINEAR THEORY OF FRACTURE, NEAR-TIP FIELDS AND INTRINSIC 

LENGTH SCALES 

The equations of motion for simple cracks, which were discussed above, are based on energy balance. We 
would like now to consider the actual helds that mediate the transport of energy to the crack tip region and 
may control the fracture process itself. LEFM predicts that sufficiently close to the crack tip the stress held 
features a universal square-root singularity, cr ^ 1/y^, where the pre-factor incorporates a universal tensorial 
function of the angle 0 and the propagation velocity v (r and 0 are polar coordinates co-moving with the tip), 
and the scalar multipliers — the stress intensity factors. Due to linearity, the displacement gradients feature 
the same singularity where is the linear elastic displacement held. 

The LEFM helds possess certain peculiar features that are worthwhile noting. For example, in mode I 
(tensile) fracture, the crack parallel component of the stress tensor ahead of the tip, cra,a,(r, ^ = 0, u), is larger 
than the tensile component (jyy{r^0 = 0,u) for every hnite crack velocity v > 0 [29143T] . In fact, the ratio 

diverges as v^cr. This feature implies that it is not easy to explain why mode 
I cracks propagate perpendicularly to the tensile loading {y) axis based on the LEFM universal helds. While 
this feature has been recognized in the past, it has not received much attention. A corollary of this is that 
there always exists a crack propagation velocity vq = ^ — q) (cf. Eq. (49) in [T] and note that Cg 

is the shear wave-speed) above which dyU^y\r^0 = ^^v) becomes negative. That is, the universal LEFM helds 
predict that tensile cracks, which ultimately involve extending — and eventually breaking — bonds ahead of 
their tips, propagate with a predominantly compressive displacement gradient ahead of their tip [32] . 

These features suggest that the LEFM universal helds are insufficient to fully explain the fracture process, 
which is supported by explicit measurements of the strain helds surrounding dynamic cracks [33] • This by no 
means implies that LEFM is wrong, this theory accurately predicts the hux of energy into the tip region and 
the universal LEFM helds pose the proper boundary conditions on any inner problem which is dehned inside 
the process zone. It does imply, though, that the length seale at which these helds are valid is quite distinct 
from the scale where fracture actually takes place. LEFM, however, is a scale-free theory that provides no 
specihc information about where it is valid and on what length scale it breaks down. This issue of the absence 
of an intrinsic length scale may have serious consequences in the context of dynamic instabilities, as will be 
discussed below. 

Recently, these issues have been addressed in the framework of the “Weakly nonlinear theory of fracture” 
mEnsii. This theory is based on a general expansion of the strain energy functional U in the form 

UiE) - ^ (trEf + II trE^ + /3i {trEf + pitrE trE^ + p^trE^ + Oi.E'^) , 


(3) 
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where E = |[(Vn) + (Vn)^ + (Vn)^(Vn)] is the nonlinear Green-Lagrange strain tensor [38] and u is the 
displacement field. Here A and /r are the standard Lame (second order) constants and {/3i,^^ 2 ,/^s} are the 
third order elastic constants [39] , which are basic physical quantities that are related to the leading anharmonic 
contributions to the interatomic interaction potential. The third order elastic constants are not free parameters, 
but rather quantities that are either calculated from a fully nonlinear elastic energy functional, if known, or 
measured directly in experiments [T]. 

As is evident from Eq. the weakly nonlinear theory takes into account the leading elastic nonlinearities 
when LEFM breaks down. The basic idea is that in the presence of the intense displacement gradients near 
crack tips in brittle materials, linear elasticity first gives way to nonlinear elasticity, rather than to dissipation. 
As such, it is an asymptotic (in space, where LEFM breaks down near the tip of a crack) and perturbative (in 
the magnitude of displacement gradients) theory based on the expansion Vn = + 0(|Vn^^^|), 

where the superscripts denote different orders in the displacement gradients, simply corresponds to LEFM 
and satisfies a new equation that depends both on the nonlinearity and The solution to this equation 
takes the form (see details in [T]) 


.( 2 ) 




[\ogr + fy^{e,v)] , u 


( 2 ) 






(4) 


and is inherently mode I (tensile) for corresponding to both mode I and mode 11. That is, the explicitly 
calculable functions /^’y^(6^,'^) satisfy {—O^v) and Note that r 

inside logr in Eq. 0 should be made non-dimensional, but this simply adds a constant and is not included 
here. 

The theory, as is clearly observed in Eq. 0, predicts a stronger than l/^/r divergence of the displacement 
gradient tensor field, and a logarithmic stretching of the crack tip opening displacement, which is 

strictly parabolic in LEFM [3]. This 1/r displacement gradient singularity has a fundamentally different status 
compared to apparently similar singularities in LEFM gnun]. In particular, its contribution to the J-integral 
vanishes identically. This, and other properties, are discussed in detail in [T] [35]. These predictions were 
quantitatively verified by direct experiments on rapid cracks, as shown in Fig. [^ In particular, it is explicitly 
shown (panel b) that for a propagation velocity for which LEFM predicts a negative displacement gradient 
dyU^\r^0 = 0,v) <0 near the tip (i.e. v>vo, where in this case 'T’o = 0-73cs), both the experimentally observed 
and theoretically predicted dyUy{r,0 = 0^v) are positive, the latter directly due to the nonlinear contribution 

dyuf\r,0 = t)^v). 

To sum up, the weakly nonlinear theory of fracture — and its experimental verification — show that elastic 
nonlinearities give rise to a stronger near tip singularity compared to LEFM and to an extensional displacement 
gradient ahead of the tip, dyUy{r^0 = {)^v) >0, even at propagation velocities where LEFM predicts otherwise 
(i.e. for v>vq^ cf. Fig. [^). The latter at least recovers the basic physical intuition that mode I fracture occurs 
under extensional deformation ahead of the propagating tip. As was mentioned above, when nonlinearities are 
taken into account for mode II cracks [42], a symmetry-breaking effect in which LEFM mode II (shear) cracks 
feature mode I (tensile) properties naturally emerges in this theoretical framework [43] . 

As the weakly nonlinear theory features two different singularities associated with two different constitutive 
relations, it predicts the existence of a new intrinsic (i.e. non-geometric) length scale, inU which characterizes 
the scale at which the two singular solutions are comparable in magnitude, i.e. roughly when | | ^ | |. 

This intrinsic length scale is a dynamic quantity (i.e. it exhibits a non-trivial velocity dependence) that 
depends on the third order elastic constants. Its length dimension is inherited from T/E, i.e. ^^zocT/E, with 
a nontrivial and not necessarily order unity pre-factor. Finally, we note that very recently an Extended Finite 
Element Method (XFEM) analysis of nonlinear fracture problems based on the weakly nonlinear theory has 
been performed mi. It was shown that using the analytic properties of the solution in Eq. 0 to construct 
crack tip enrichment functions provides an accurate and efficient tool for solving finite strain fracture problems. 

The weakly nonlinear theory of fracture, which invokes three directly measurable third order elastic constants, 
provides a general and systematic way to go beyond LEFM at smaller scales near the tip. It naturally resolves 
apparently paradoxical physical features of the LEFM universal fields and gives rise to a new intrinsic length 
scale, which is missing in LEFM. As will be discussed soon, the latter plays an important role in dynamic 
instabilities. 
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FIG. 3: (a) The measured yy-component of the displacement gradient tensor eyy{r,0 = 0) = dyUy{r,0 = 0) ahead of a 
crack propagating at i’ = 0.20cs (circles) compared to the LEFM prediction (dashed line) and the weakly nonlinear theory 
prediction (LEFM supplemented with Eq. @) (b) The same for r' = 0.78cs. Note, in particular, that for v>vo = 0.73cs, 
LEEM predicts Syy{r, 0 = 0)<0 sufficiently close to the tip, while both the experiment and the weakly nonlinear theory 
show that eyy{r,0 = 0) >0, as is physically expected in mode I fracture, (c) The measured crack tip prohle (0^(r, 0 = ±7r) 
vs. (I)x{r,0 = 7r) for r' = 0.53cs, where (p{x)=x + u{x)) (circles) compared to the best ht of the parabolic LEEM prohle 
(dashed line) and that predicted by the weakly nonlinear theory (solid line). The effect of the logarithmic stretching 
term in Eq. @ is evident as the tip is approached. Adapted from [34], where additional details can be found. 


IV. CRACK INSTABILITIES 
A. The micro-branching instability 

Historically, early measurements showed that a running crack did not seem to be described by the equations 
of motion described in Sect.jllj In particular, these measurements suggested that even as robust a prediction as 
the limiting velocity of a dynamic crack seemed to be off; cracks in brittle materials rarely surpassed even 60% 
of cr. These observations suggested that something was seriously wrong with our fundamental understanding 
of rapid crack dynamics, but it was unclear where the root of the problem lay. One could always “blame” 
the velocity dependence of the fracture energy, T{v), but then one had to explain why, in brittle amorphous 
materials T{v) needed to increase precipitously in order to prevent Eq. 0 from allowing a crack to approach 
Cr. Even more puzzling was the fact that this behavior occurs in a variety of materials whose molecular 
dissipative mechanisms are expected to be wholly unrelated; these materials ranged from brittle (cross-linked 
and uncross-linked) polymers to glasses. 

The answer to the apparent paradox suggested above is that at high velocities, the crack is no longer a 
“simple crack”. Instead, as first observed in mills], a crack’s structure at relatively high propagation velocities 
is entirely different from that envisioned in the picture of simple crack; high-speed photographs of rapid cracks 
showed that the structure of the crack tip was anything but simple. Instead of the propagation of a single simple 
crack, these experiments showed that rapid cracks are composed of an ensemble of simultaneously propagating 
cracks. Later experiments [46], suggested that the break-up of a simple crack to multiple cracks occurred at a 
critical velocity, at which a single crack became unstable. This instability, which became known at the “micro¬ 
branching” instability, was later shown to be characteristic of many brittle amorphous materials miEIlllSHSSI- 
A comprehensive review of some of the early work on the micro-branching instability can be found in m- 

What is the micro-branching instability? For low crack velocities, the crack propagates smoothly and the 
fracture surfaces formed are smooth and mirror-like, with the only apparent structure caused by material defects 
interacting with the singular front [54l|56]. At a critical velocity of Vmb — ^-^CR, the crack undergoes a dynamic 
instability. As shown in Fig. rapid fluctuations appear in the instantaneous crack velocity and, at the 
same time, non-trivial structure emerges on the fracture surface. The fracture surface structure is accompanied 
by frustrated crack-branching, as multiple microscopic crack branches (“micro-branches”) are observed below 
the fracture surface. As the crack velocity increases beyond Vmb^ both the length and density of the micro- 
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FIG. 4: A comparison between the instability in both soft polyacrylamide gels (left) and soda-lime glass (right). Once 
a crack attains a critical velocity of i’ ~ 0.4ci?, a single “simple” crack may become unstable to the micro-branching 
instability. At this point (top) the instantaneous velocity of the crack undergoes violent oscillations that correspond 
to the formation of subsurface frustrated micro-branches (center), where side views {xy plane) of micro-cracks in both 
materials are presented [x is the propagation direction, y is the loading direction and z is the direction along the 
crack front). Micro-branches are very similar in both materials, characterized by a power-law functional form, (bottom) 
Photographs of the resulting fracture surface [xz plane) formed by cracks propagating from left to right at u~0.5ci^. The 
chains of structures on the fracture surface correspond to chains of micro-branches that are aligned in the propagation 
direction x and highly localized in the z direction. The vertical dimensions of the photographs are about 0.5mm. 
The cartoon in the center describes the geometry; the xz fracture surface is formed by a propagating cracks whereas 
micro-branches (dotted lines) extend below the surfaceAdapted from [T]. 



FIG. 5: A series of photographs of the quasi-parabolic prohle of running cracks in polyacrylamide gels. The two lower 
left panels depict the parabolic crack tip opening prohle of “simple” cracks propagating just prior to the onset of the 
micro-branching instability (corresponding to the green dot on the velocity vs. crack length hgure in the upper right 
corner). Right panels: 4 prohles of cracks undergoing the micro-branching instability, corresponding to the red dots in 
the velocity measurements. Note the small micro-cracks formed at the tip of the main crack. Adapted from [T]. 
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branches increase with v. The increased length and density of the micro-branches gives rise to correspondingly 
rougher fracture surface features, larger velocity fluctuations, and a sharp increase of the fracture energy that is 
proportional to the total micro-branch length [49j[57]. As Fig. [^demonstrates, these characteristic features are 
independent of the brittle material; shown are two extremely different classes of material, poly-acrylamide gel 
(an aqueous elastomer) and soda-lime glass which exhibit identical characteristic behavior [53]. Note the use 
of brittle gels has enabled experiments to probe the fracture process in unprecedented detail, as demonstrated 
in Fig. 

Over the past two decades, there have been a number of notable attempts to describe the micro-branching 
instability by supplementing or extending LEFM in various ways. These include phase-field models [58f[64] . 
phase transformation models [65|, cohesive-zone models misniis], lattice models IS7H73], models based on 
the “Principle of Local Symmetry” [741476] , energetic bounds on crack branching [TMS] and models based on 
non-linear constitutive behavior near the crack tip [801482] . Although many of these models exhibited some 
properties that were qualitatively consistent with the experimental observations, they are inherently deficient 
because they focussed on 2D (i.e. they actuaiiy studied macro-branching, which appears as a tip-spiitting 
instabiiity). As wiii be shown beiow, micro-branching appears to be an intrinsicaiiy 3D instabiiity that cannot 
be satisfactoriiy described as a 2D phenomenon. 

On the conceptuai ievei, we note that most of previous works have considered the near-tip region as a passive 
region that both reguiarizes the singuiar fieids predicted by LEFM and accounts for the dissipative processes at 
the tip, but otherwise piays no other active roie in the dynamics. As wiii be shown next, recent deveiopments 
may suggest that understanding crack instabiiities may entaii the introduction of new physicai ingredients (e.g. 
iength scaies) in which the near-tip region plays a more active role. 


B. The oscillatory instability 

By manipulating the geometry and thickness of the brittle gels mentioned above, it is possible to statistically 
suppress the micro-branching instability in sufficiently thin samples (200 — 300/im) [83] [84]; while micro¬ 
branching can occur for all velocities 0.3cs < v < O.Ocg, the probability of exciting them decreases for both 
thin samples and large crack accelerations. When micro-branching is suppressed, simple brittle cracks can 
reach extreme velocities. In such cases single crack dynamics, as described by LEFM (Fig.[^, can be obtained 
for velocities approaching cr. 

Can these cracks actually reach their asymptotic velocity? It turns out that this is prevented by a new 
instability. Experiments in gels have found that at velocities surpassing about 90% of the shear wave-speed, c^, 
a straight crack becomes unstable to sinusoidal path oscillations [83] . This instability is particularly intriguing 
since it involves a finite instability wavelength at onset. The value of the observed wavelength was not dependent 
on either the system geometry or the loading conditions. This suggests the existence of an intrinsic scale that 
can not exist in linear elastic solutions for cracks, which are scale-free. 

This oscillatory instability should be distinguished from apparently similar oscillatory cracks observed in 
rubber [85]. The oscillatory cracks in rubber were intersonic {cg <v<Cd) [86], entailed bi-axial loading, and 
were subjected to very large background strains (which make the material response nonlinear everywhere). 
The oscillatory cracks in the elastomer gels were subsonic and were observed under uniaxial loading of such 
magnitude (^ 10%) that the deformation away from the tip was predominantly linear elastic. 


C. Understanding 2D instabilities 


In Sects, [n] and |III| we discussed simple tensile cracks that propagate along straight paths in effectively 
2D bodies. By so doing, we excluded spontaneous symmetry-breaking instabilities from the discussion. Such 
dynamic instabilities, however, are rather prevalent and pose some of most intriguing open questions in the 
field of dynamic fracture. In sections IV A| and |IVB we described two different instabilities that take place in 
fast fracture. 

The oscillatory instability described in Eig. [^is, in a sense, the simpler of the two instabilities. Once excited, 
the oscillatory instability breaks the translational invariance of a crack in the propagation (x) direction. The 
fracture surface created during the crack’s motion is wavy but mirror-like; there is no apparent structure along 
the fracture surface. Thus, we can consider the crack front to be translationally invariant along the direction 
transverse {z direction) to the crack’s propagation direction. It is, therefore, sufficient to idealize the crack tip 
as an oscillating point within an effectively 2D medium. 
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FIG. 6: (a) Typical photographs of the xy profiles of fracture surfaces at the onset of the oscillatory instability; from top 
to bottom: gels with shear moduli /i = 36kPa, 143kPa, 168kPa. (b) The oscillation wavelength, Aosc changes significantly 
with the material, as characterized by y. Symbol colors as in (b). Adapted from [84] . 


In contrast to the oscillatory instability, when the micro-branching instability is excited, distinct structures 
are generated along the fracture surface that are highly localized along both the x and 2 : directions (see, for 
example, the lower panels of Fig. [^. As the crack front is not translationally invariant along the ^ direction, 
it may not be sufficient to consider a crack propagating within a 2D medium to understand micro-branching. 

We now consider the oscillatory instability in a 2D medium. As mentioned above, the most interesting 
physical property of this instability, beyond its existence, is that its wavelength does not scale with the geometric 
dimensions of the sample. This implies that it cannot be explained by LEFM, as it possesses no intrinsic length 
scale beyond geometric length scales. What is then the origin of this instability and the emerging intrinsic 
length scale? 

To address this question, and dynamic fracture instabilities in a broader context, we need to consider first 
a basic aspect of the current theoretical status of the field. When cracks propagate along straight paths 
(or other predetermined paths), the quantity of interest is the rate of crack growth, which involves a scalar 
relation based on energy balance. In the most general case, however, the path is an intrinsic property of 
the dynamic fracture problem and is selected self-consistently as in a large class of moving free-boundary 
problems (see, for example, [87H92] b Consequently, we need an equation of motion for crack tips (in 2D) 
and crack fronts (in 3D) that will allow the determination of both the direction of propagation and its rate. 
Such equations are intrinsically non-scalar in nature and involve physics that go significantly beyond energy 
balance alone. In particular, such equations are essential in the context of instabilities as without them one 
cannot even mathematically pose the question of stability, which amounts to perturbing solutions of existing 
and well-defined differential equations. We currently lack well-established crack tip/front equations of motion 
for general, non-straight, dynamics. 

In 2D and under quasi-static conditions, v<^cr, it has been shown that the Principle of Local Symmetry 
[93] offers good approximations for crack paths |94Hl()()j . This approximation suggests that under quasi-static 
conditions, and under the assumption that no intrinsic length scales exist near the tip, crack propagation 
along smooth paths will annihilate any shear (mode II) component of the fields ahead of the tip. Hence, this 
path selection criterion proposes that conditions of local tensile symmetry are almost maintained at a crack’s 
tip, which mathematically amounts to KnjKi ^1. As the stress intensity factors are generically non-local 
functionals of the path, the latter is in fact an integro-differential equation. At present, to the best of our 
knowledge, there is no evidence that the Principle of Local Symmetry offers good approximations also under 
fully dynamic conditions. In fact, some recent experiments might indicate otherwise uni]. 

Recently, a dynamic crack tip equation of motion, which takes into account the existence of an intrinsic 
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FIG. 7: (a) The scaled critical velocity for the onset of the oscillatory instability, I’c/cg ~0.9, is constant. Vc is defined as 
the maximal velocity prior to the instability onset in each material. Symbol colors correspond to the legend in (b). (b) 
Comparison between an experimental measure of the length scale associated with elastic nonlinearities 5{vc) (measured 
at the critical velocity I’c, see [84] for details) and the oscillation wavelength Xosc- Note that the different combinations 
of 11 and r are used to produce these 15 independent measurements. The dashed line is a guide to the eye. Adapted 
from [84] . 


length scale £ni and the possible high propagation velocities, was proposed [82]. The basic idea is that in the 
presence of a finite length scale near the tip, causality implies a time delay between the singular fields 
driving the crack (as described by LEFM) and the actual crack tip. Taking into account this time delay along 
with scaling arguments and symmetry considerations [102] . a dynamic crack tip equation of motion has been 
constructed [82] . 

This equation was used to study the linear stability of straight crack propagation against sinusoidal pertur¬ 
bations, with the help of the Willis-Movchan perturbation theory for the stress intensity factors [10311105] . This 
analysis predicted an oscillatory instability of mode I cracks to infinitesimal out-of-plane perturbations above a 
high propagation velocity of 'T’c —O-Scg, with the following properties: (i) The normalized critical velocity Vc/cs 
is nearly universal (it exhibits a weak dependence on Poisson’s ratio u) (ii) The wavelength of oscillation is 
proportional to the intrinsic length scale /E [84]. These predictions were recently tested by extensive 

experiments in which F and E were varied over a substantial range [84]. As Fig. shows, a favorable agreement 
between theory and experiments was demonstrated. 

We believe that these results are important because they show that the process zone is not a passive region 
responsible only for nonlinearity and energy dissipation, but rather an active region whose accompanying length 
and time scales can significantly affect macroscopic fracture dynamics, in particular crack instabilities. The 
importance of a length scale proportional to T/E has also been highlighted in experiments on the quasi-static 
fracture of hydrogels, where it was shown to play a decisive role in a mode I cross-hatching instability [ins], 
in a wetting-induced branching instability and in a mixed-mode I+III (mode III corresponds to anti¬ 

plane shear fracture) echelon instability [108] . We suspect that the incorporation of intrinsic length scales into 
fracture theory may be essential to developing equations of motion for crack tips/fronts and to understanding 
dynamic instabilities. Whether or not these ideas are also relevant to the crack oscillations observed in rubber 
[HiiHaiing], where the applied strains are extremely large and intersonic propagation velocities were reached, 
remains an open question. 


V. THREE-DIMENSIONAL TENSILE CRACKS: MICRO-BRANCHING, FRONT WAVES AND 

CORRUGATIONS 

In the previous section we discussed a way to understand the oscillatory instability based on a dynamic crack 
tip equation of motion that takes into account the length scale proposed by the weakly nonlinear theory of 
fracture mechanics. It was shown that a crack propagating in a 2D medium will lose its stability at extremely 
high velocities (^ O.Qcg) to infinitesimal out-of-plane perturbations. We are still faced, however, with the 
challenge of explaining the micro-branching instability. Micro-branching, on the face of it, appears to be a 
much more difficult problem. Once excited, this instability appears to be intrinsically 3D, as Fig. [^suggests. It 
had been discovered in earlier work [53] that the micro-branching threshold, Vmb^ could vary above a minimum 
value of about — The observed values of Vmb were shown to statistically depend on the value of the 

acceleration prior to the instability initiation. In addition, the transition to micro-branching was observed to 
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FIG. 8: Phase diagram of the micro-branching instability. The micro-branching instability is triggered when a small 
Kii component is present. For each propagation velocity there exists a threshold value, Kf}(v). When Kn > Kf}{y) 
micro-branching is observed. Shown are values |iFjJ/iFj |,which are the values obtained for each v immediately prior 
to micro-branching. The critical boundary linearly extrapolates (dashed line) to zero at i’ ~ 0.87cs, the value of the 
measured threshold velocity for the onset of the oscillatory instability [84]. Colors denote different experiments. Adapted 
from [lOlj . 


be hysteretic in v. 

These observations suggested that the micro-branching instability involves an activated process that could 
be triggered by noise above the minimum value of If the acceleration is low, micro-branching is always 
observed near the minimal value of Vmh^ as eventually the crack front will encounter a noise source that is 
sufficiently large to trigger the instability. In very thin gel samples, where “sources of noise” along the z 
direction are geometrically reduced, and at high crack accelerations, micro-branching was seen to occur in a 
relatively large range of velocities, 0.35cs <0.9cs. This is the reason that the oscillatory instability could 

be, at all, observed prior to the onset of micro-branching [831(84] . 

Very recent experiments, however, suggest that the micro-branching instability may be closely related to the 
oscillatory instability m\- These experiments directly measured the 2D strain fields surrounding the tips of 
rapidly moving cracks in thin gels, as these cracks approached the instability over a wide range of Vmb values. 
Prior to the onset of micro-branching, the cracks were “simple” so that the media could be considered as 2D. It 
was discovered that the strain fields of all of these simple cracks prior to the onset of the instability contained 
a fluctuating shear component that could be characterized by the textbook Ku fields [3] . Examining the value 
of Kjj immediately prior to the instability onset, these studies showed that there exists a velocity-dependent 
threshold value, for the onset of instability. 

iFjj as a function of v is presented in Fig. The data, when extrapolated to Kfj{v) 0, yield a zero 
threshold propagation velocity of O.QCg, which is precisely the critical velocity for the onset of the oscillatory 
instability. Guided by local symmetry arguments [93], we would expect any finite Ku value to cause a crack 
to curve away from the initial symmetry axis imposed by the loading. The data therefore suggest that two 
distinct “phases” can coexist in a hysteretic region 0.3cs <'r’<0.9cs, where both straight cracks and cracks that 
locally curve out of the symmetry plane can be simultaneously stable. Local crack curvature will only occur 
upon a sufficiently large (finite) local shear perturbation, as shown in Fig. Moreover, in this bi-stable region, 
both straight and curved cracks can coexist within different sections along the crack front. This bi-stability 
provides an explanation for the localized nature of the micro-branches in the ^-direction. The bi-stability that 
enables this localization disappears at the critical velocity 0.9cs, where \Ku\^^ = 0. At this point the entire 
crack front must coherently curve out of plane and the oscillatory instability ensues. 

Insofar as propagation along the symmetry axis is concerned, it is clear that straight sections of the crack 
front will outrun the curved sections. As a result, we expect that the crack front will eventually revert to a 
single fracture plane; pinching off the curved section. This “frustrated curved section” is what we identify as 
a micro-branch beneath the fracture surface. Such a pinch-off scenario been recently observed experimentally 
by direct visualization of crack fronts in the xz plane as micro-branching takes place [llOj . The micro-branch 
locally delays the crack front causing local crack front curvature within the xz (fracture) plane. This initial 
curvature was seen to continually deepen, eventually developing into a cusp-like shape. Cusp formation within 
the fracture plane was coincident with the “death” of the micro-branch. 
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Micro-branching is, therefore, an intrinsically 3D instability that is intimately related to the crack front 
dynamics, and hence defies any attempt to understand it in a purely 2D framework. There are several other 
basic phenomena that necessitate a 3D framework of fracture dynamics. In the quasi-static regime, simulations 
coupled with analytic work |111I4115] has shed light on a crack front instability under mixed-mode I-j-III 
conditions.Recently, helical perturbations have been shown, in 3D simulations, to develop into segmented crack 
fronts |11614118j . Related segmented crack front structures have been observed in mode I [10611119j and mixed¬ 
mode I-f III fracture of of both “standard” |12Ql[125j and soft materials |lQ8j . 

In the dynamic regime, crack fronts were shown to support propagative modes along them. The most well- 
established of those are in-plane crack front waves, which were predicted through perturbative mode I numerical 
[231 mi and analytic [22] calculations in the late 1990’s. In this case, the crack front is assumed to remain 
within its tensile symmetry plane, but to experience in-plane shape perturbations. In the framework of LEFM 
and a linear perturbation theory [22ll24] . perturbations to the crack front were shown to propagate along the 
front without attenuation (as long as T is independent of v) at nearly cr^ for any v. 

Subsequent experiments on various amorphous brittle materials [25H28] indeed revealed the existence of 
localized modes that propagated along the crack front at nearly cr^ when triggered by an interaction with 
localized material inhomogeneities (asperities) or by spontaneously generated micro-branches. The experimen¬ 
tally observed propagative modes were different from the theoretically predicted in-plane front waves in two 
important aspects: (i) They featured an out-of-plane component, which left marks on the fracture surfaces 
and, in fact, enabled post-mortem measurements of crack front waves (in-plane perturbations do not leave any 
post-mortem fractographic signature) (ii) They were clearly nonlinear, exhibiting soliton-like properties [25] . 
While similar fracture surface markings can be generated by the interaction of shear waves with moving crack 
fronts |12614128] , recent measurements in both brittle gels [53[ IllQj and the rapid fracture of rock [28[ 1129] have 
strengthened the evidence for the existence of front waves in dynamic fracture. 

Treating out-of-plane front perturbations is more difficult than in-plane ones because the front interacts 
with out-of-plane perturbations which are locked in the front’s wake and because energy balance alone is 
insufficient to determine the front dynamics (information about the direction of propagation is essential as well, 
see discussion above on crack front equations of motion). Only very recently, based on the Willis-Movchan 3D 
LEFM linear perturbation formalism nni and the Principle of Local Symmetry, has the out-of-plane stability 
of crack fronts been considered [13011131] . In particular, a minimal scenario in which linearly unstable out-of- 
plane crack front corrugations might emerge (above a critical front propagation velocity), has been discussed. 
Within this scenario, a critical velocity has been calculated as a function of Poisson’s ratio and corrugations 
have been shown to propagate along the crack front at nearly the Rayleigh wave-speed [m]. 

To linear order in the deviation from a straight crack front, in- and out-of-plane perturbations are decoupled. 
Taken together with the nonlinear nature of the experimentally observed front waves, nonlinear coupling be¬ 
tween in- and out-of-plane front dynamics appears to be an essential next step. Due to the technical difficulties 
involved, numerical approaches such as phase-field models may be useful in achieving this In addition, 

the effect of going beyond the Principle of Local Symmetry, e.g. in the spirit of [82], should be systematically 
explored. Possible connections to the much studied problem of fracture surfaces roughness [132] . where dynamic 
effects are traditionally neglected urn, should also be investigated. 


VI. THE DYNAMICS OF FRICTIONAL INTERFACES: ON THE RELATION BETWEEN 

FRICTION AND FRACTURE, AND MORE 

In this chapter we describe recent work that suggests that our understanding of many aspects of “common¬ 
place” friction may need to be re-examined. This work suggests that dynamic fracture processes play a major 
role in frictional motion - even when these processes appear to be extremely slow. The key to understanding 
how dynamic fracture comes into play is to consider the spatially extended interface that separates two bodies 
in frictional contact. In dry friction, this interface is composed of a large ensemble of rough and interlocking 
discrete contacts, whose total area of contact is orders of magnitude smaller than the nominal contact area of 
the contacting bodies [13411135j . In order for frictional motion to initiate, these contacts must be “broken” or 
rupture. As we know, fracture of a spatially extended entity is not simultaneous, but is generally mediated 

by crack-like fronts. Below, we describe observations of such rupture fronts, which will be shown to be closely 

related to classic dynamic shear (mode II) cracks. It is intuitively obvious that these rapid rupture modes are 
also closely related to earthquakes, which mediate the relative motion of frictionally coupled tectonic plates. 

The failure of macroscopic frictional interfaces has for centuries been thought to be governed solely by 

“friction laws”, with the classic example of Amontons-Coulomb law of friction Fg = where Fg is the 

friction (shear) force and F^ is the normal force. The proportionality constant is either the static friction 
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FIG. 9: Different types of rupture fronts in friction. Measurements of the real contact area reveal three qualitatively 
different types of fronts that “fracture” the contacts defining the interface that separates two sliding bodies. Shown 
is the space-time variation of the real area of contact along a 200mm long frictional interface. The sharp reduction in 
the contact area occurs at the tip of the propagating fronts. Shown are typical examples of (top) slow (center) fast 
(sub-Rayleigh) and (bottom) supershear (i; > Cs) fronts, (adapted from [141] b 


coefficient fig or the dynamic one both are considered to be characteristic properties of the materials 
forming the frictional interface. This class of friction laws has been generalized to spatially varying and more 
physically faithful constitutive laws [13614140] . The latter, termed rate-and-state friction laws, incorporate the 
effects of additional dynamic variables such as the local sliding velocity (“rate”) and the local average contact 
lifetime (“state”) on friction. 

A new generation of experiments over the past decade has revealed that rapid rupture fronts drive frictional 
motion. In other words, the transition from “static” to “dynamic” friction has been shown to be mediated by 
different types of crack-like fronts. These fronts are the vehicles that trigger the onset of macroscopic frictional 
motion, bridging the gap between the microscopic interactions that define local frictional resist^ce and the 
resulting macroscopic motion imbued in the slip of large bodies [14244145] . As shown in Fig. these new 
experiments have revealed various rupture modes that are loosely classified as follows: slow ruptures, which 
propagate far below the material wave speeds [m ESHllSl, “sub-Rayleigh” ruptures [Mil 1^ UlTl [T49] 
that propagate up to the Rayleigh wave-speed cr^ and “super-shear” rupture modes that surpass the shear 
wave-speed Cg [14111146] I149H152] . It is important to note that analogs of all of these rupture modes have been 
documented in natural earthquakes [144111451115311156] . 

Are the rupture dynamics of an interface related to the values of the global “friction coefficients” that 
characterise the interface strength? Recent experiments using PMMA and granite [152] have measured 
the static friction coefficient, i.e. the ratio between the applied shear and normal forces at the onset 
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FIG. 10: Count of events with varying /Xs, from experiments conducted under widely varying loading conditions. These 
range from pure edge loading to near uniform application of the shear force, Fs, to the system. Note that this plot is 
not a histogram and solely reflects the number of experiments performed whose loading conditions yielded fig for each 
range. The color bars (bar width = two standard deviations) in the plot depict sub-populations of events from single 
experiments, each with different (specific) loading conditions. Adapted from |157j . 


macroscopic sliding. These studies found that sliding can be initiated over a wide range of static friction 
values, indicating that jUg is certainly not solely a characteristic material property but is, instead, dependent 
on how a given system is loaded Study of the rupture dynamics leading to the broad range of friction 

coefficients, in fact, revealed that the loading conditions of the systems and the subsequent rupture dynamics 
were closely correlated to the value of jUg |157j . These experiments found that fig was as well defined as the 
loading conditions. As long as the way in which forces were applied to two blocks in frictional contact was 
kept constant, fig was extremely well-defined. Once the loading configuration was changed, the measured value 
of fig for the same two blocks could be varied significantly and controllably. This surprising result has been 
confirmed in recent theoretical and numerical studies |158H161] . 

In addition to analytical work, recent numerical work [1581 nsni usa css] together with simple models 
describing the contact dynamics of ensembles of discrete contacts have been successful in describing many 
aspects of frictional cracks |148lll64tll67j . These results have provided important new insights into the influence 
of both the discrete nature of the contacts that compose the frictional interface and the influence of different 
“friction laws” that couple the elastic medium to the interface. 

How close to “ordinary” fracture might the rupture processes shown in Fig. I^be? Let us consider two elastic 
bodies in frictional contact, subjected to remotely applied shear stresses. While in an intact homogeneous 
system the Principle of Local Symmetry would predict that any crack should rotate so as to annihilate any 
shear (mode II) component near its tip |168j . the presence of a weak friction plane (compared to an intact 
bulk) allows rupture fronts to propagate under predominantly shear conditions along the weak plane. Are 
these rupture fronts related to classical mode II cracks? The surfaces left behind classical mode II cracks are 
assumed to be traction-free. While the propagation of rupture fronts along frictional interfaces is definitely 
accompanied by stress reduction (“stress drop”), it is clear that the frictional stress behind the propagating 
front is finite, and hence the physical situation is evidently different from the traction-free boundary conditions 
assumed in classical mode II crack solutions. 

The two problems (i.e. of frictional rupture fronts and mode II cracks), though, are closely related when the 
following two conditions are met: (i) The residual stress left behind the rupture front is constant, a situation 
commonly assumed in the context of earthquake dynamics |142[ 11431116914171j (ii) The width of the front (the 
“cohesive zone”), i.e. the length scale over which the residual stress is reached when the front passes a point 
at the interface, is small compared to other length scales in the system. In particular, under these conditions 
the singular mode II fields of LEFM remain valid (in this case the singular Kn field mediates the transition 
from the initial shear stress to cr^) [T7^[T73] . 

These predictions have been directly tested in very recent experiments on PMMA blocks, where the real-time 
strain tensor at points adjacent to the interface was measured during rupture propagation [mi. In parallel, the 
instantaneous rupture tip location and velocity were measured. These techniques, which were used successfully 
to describe the strain fields in mode I fracture HZS], provided a first-time quantitative description of the strain 
fields surrounding the tips of these propagating rupture fronts. Measurements were performed for rupture 
fronts over a broad range of propagation velocities x;, 0.04ci^ <x; <0.97ci?. 
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Typical results of these experiments are shown in Fig. El where a comparison of the different components of 
the strain tensor to the classic mode II singular solution for rupture fronts propagating at both slow velocities 
{v < 0.3ci?) and velocities approaching the Rayleigh wave-speed, is presented. The surprisingly good 
agreement demonstrates that, at least in this material, the two conditions listed above are satisfied, and 
consequently the fields near the propagating rupture fronts are quantitatively described by the singular solutions 
predicted for shear (mode II) cracks in the framework of LEFM. It is interesting and important to note that the 
agreement to the LEFM singular fields was achieved with a velocity-independent fracture energy F, suggesting 
that the well-documented rate dependence of friction in PMMA [nni does not strongly influence the effective 
fracture energy. 
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FIG. 11: (A) Measurements of the strain tensor, at a fixed spatial location above the frictional interface (see panel 
C for details) during the passage of a rupture front. The initial strain was subtracted from Sxx and Syy and residual 
strain from Sxy Colors represent different, relatively small, front velocities {0.05cr<v <0.24cr). Corresponding LEFM 
predictions for plain strain boundary conditions are plotted in black {v = 0.2cr). The fracture energy F~l.lJ/m^ is the 
sole free parameter. (B) Sij for u~0.96ci?. The LEFM solution (black), plotted as in panel A, still describes the larger 
amplitudes and strong oscillations of rapidly propagating ruptures with the same fracture energy as before. (C) Shear 
strain variations, Exy , surrounding the rupture tip predicted by LEFM for v = 0.96ci?. The blue square denotes the 
strain gage rosette measurement location relative to the approaching rupture tip. Note the strong angular dependence 
that drives the violent oscillations in panel B. Adapted from mi¬ 
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The experiments described in Fig 

with mode II cracks. Let us now consider the question of rupture nucleation. 


suggest that propagating sub-Rayleigh ruptures can be identified 

This question is of great 


importance both in general and in the particular context of earthquake nucleation (see, for example, [17611180] . 
If nucleation takes place as a Griffith-like process analogous to conventional crack initiation [181] , the nucleation 
length Lc — the critical size above which a slowly creeping patch becomes unstable — is determined by the 
physical conditions near the rupture tip and depends on the pre-stress level. Alternative friction-law dependent 
nucleation scenarios, however, are possible. For example, rate-and-state friction models [na Ei uTTHiHoi 
associate Lc with a spatiotemporal friction instability of quasi-statically extending “creep patches” that is 
independent of the pre-stress level [TT^ITM] . It is currently unclear which nucleation scenario, if any, is more 
physically realistic and, if so, under what conditions. 

What is the origin and nature of the experimentally observed “slow fronts” fTTniTM] ? Recent theoretical 
work has suggested that they may be a direct result of the qualitative form of the friction law at the interface; 
the frictional constitutive law can also affect the existence and nature of propagating rupture fronts. To see this. 
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FIG. 12: A steady sliding friction curve fi(Au), where Aii is the sliding velocity, which is typical of many materials [183] . 
The curve exhibits a non-monotonic behavior with a minimum after a small sliding velocity maximum. The three 
horizontal lines correspond to different levels of applied shear stress. A stable homogeneous solution is marked by a red 
circle and an unstable one by a green square. 


consider a frictional interface separating two infinitely long elastic blocks of height H. When the interface is 
described by a conventional velocity-weakening friction law, e.g. a friction law in which the frictional resistance 
decreases from to fid with increasing slip velocity Aff, steady-state propagating rupture front solutions can 
be found under some fixed shear displacement loading conditions. These solutions are analogous to those found 
for tensile cracks in a strip [3], where the only difference is the role played by the residual stress cr^. On the 
other hand, no steady-state solutions can be found under shear stress controlled conditions; as in the analogous 
Mode I problem of constant stress, frictional cracks should continually accelerate to cr (as seen in Fig. [^. 

What happens for a non-monotonic friction law? Such a friction law, where the steady-sliding frictional 
resistance is a non-monotonic function of the sliding velocity, with a maximum at very small velocities and 
a higher velocity minimum, is depicted in Fig. This form is representative of a broad range of materials 
m- Consider again two infinitely long elastic blocks, now separated by an interface described by a steady 
sliding frictional resistance as in Fig.[^ A steady state propagating front is an object that smoothly connects 
two stable solutions in a homogeneous system, once the system is perturbed (e.g. via a stress gradient) 
and one of these solutions loses stability to the other. In the context of Fig. stable solutions correspond to 
velocity-strengthening branches (i.e. when fi increases with increasing Ait) and unstable (accelerating) solutions 
correspond to velocity-weakening branches (i.e. when fi decreases with increasing Ait). 

When the applied shear stress corresponds to the lower horizontal line in Fig. below the minimum of 
the friction curve, there is only one stable homogeneous solution (marked by a red circle) and consequently 
no propagating front solutions exist. The situation is qualitatively different when the applied shear stress 
corresponds to the second horizontal line, which touches the friction curve at its minimum. In this case, there 
are two homogeneous solutions and a propagating front solution connecting the two is possible. Such a stably 
propagating solution was shown to exist [1651118211184j , where the characteristic slip velocity corresponds to the 
minimum of fi{Au). The solution’s propagation velocity could be significantly smaller than the elastic wave- 
speeds and independent of them [TMirrsiiiTM] . Further increase of the applied shear stress (e.g. the upper 
horizontal line) produces a continuous spectrum of propagating front solutions with increasing propagation 
velocities, which eventually saturates at the relevant wave-speed 11841 . At stresses above the minimum of the 


friction curve, an unstable homogeneous solution also exists (green square in Fig. 12). 

Such propagating front solutions do not exist for friction laws that monotonically decrease with increasing 
slip velocity. This is an example in which the properties of the friction law have qualitative implications for the 
existence of steady state rupture fronts. Transient rupture fronts, which emerge in more complicated situations, 
may be considered to be short-lived excitations of these steady state propagating rupture fronts |165[ 118211184j . 
The existence of these solutions may therefore provide a possible explanation for slow interface ruptures. 
The possible roles of non-monotonic friction laws in naturally occurring slow earthquakes have been recently 
explored |1861[19Qj . While other mechanisms for the emergence of slow rupture have also been proposed. 
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they generally involve invoking other physical ingredients such as higher dimensionality, spatial variation of 
the friction law and stress heterogeneities dilation and pore pressure effects |193j . and a discrete 

description of asperities, including stochasticity of various local quantities |167j . 

The frictional constitutive law can also give rise to crack-like objects which are not classical mode II cracks. 
For example, in the presence of viscous-friction (i.e. when the frictional resistance is linear in the sliding velocity, 
which may be relevant to both dry and lubricated frictional interfaces) the order of the displacement gradients 
singularity near the tip may differ from the common yT singularity and depend on the friction law |I94[ II95j . 
Furthermore, interesting physical properties such as a healing instability, a static threshold equivalent to the 
Griffith threshold and dissipation which is dominated by friction in a large interfacial region away from the tip 
region can emerge in certain situations [1941 II95j . Such frictional shear cracks were shown to be related to the 
experimentally observed self-healing slip pulses along an interface between a gelatin gel and a glass |I96I[I98] . 
These works established a clear connection between novel crack-like objects and frictional dynamics. 

To this point we have described some recent insights obtained in the study of both sub-Rayleigh {v < cr) 
and slow rupture fronts. We now briefly describe ruptures that propagate at velocities that are faster than 
Cr. These are commonly referred to as “super-shear” ruptures. Super-shear cracks have long been considered 
in the fracture literature 0 UMl Eon], but until they were observed in the beautiful experiments by Rosakis 
and coworkers they were often considered to be of purely theoretical value. Since this first observation 

in weakly bonded homogeneous blocks that were ballistically loaded, mode II super-shear ruptures have been 
observed numerically |2QIH2Q3] as well as in a variety of other scenarios that include frictional interfaces and 
quasi-static loading 

These recent experiments have revealed that super-shear ruptures are not rare at all, but appear in many cases 
when the initial energy stored in the surrounding material is much higher than that needed to drive a running 
crack. When loaded ballistically, super-shear states are observed only for relatively high projectile velocities 
mu, whereas in frictional sliding super-shear ruptures are initiated when the initial shear stresses are very 
large (e.g. when the “friction coefficient” is significantly greater than /j.s—0.6 in PMMA |I4Ij or granite |I52j b 
In addition to laboratory experiments and simulations, there is growing evidence that super-shear propagation 
can occur in earthquakes |I56[ I2Q4I[2Q7| which are often associated with extreme damage. 

While their existence in no longer in doubt, a more precise characterization of super-shear ruptures is needed. 
It is not yet clear, for example, if the observed super-shear ruptures along frictional interfaces indeed correspond 
to “textbook” super-shear solutions of mode II fracture (in analogy to the characterization performed in Fig. 


The nucleation of super-shear ruptures is an important and still open question. While a number of scenarios 
for the nucleation of such modes have been suggested [nanMiiinniEni], it is not yet entirely clear how the 
transition from sub-Rayleigh to super-shear ruptures comes about. A common theme among these different 
scenarios is that sufficient elastic energy becomes available in close proximity to a frictional interface. This has 
been shown to take place either via the interaction of the crack front with localized perturbations |2Q8j . via 
radiation from an existing sub-Rayleigh front or simply due to sufficient pre-existing strain 

energy surrounding the interface prior to nucleation |I4Ij . This last scenario has been suggested as a necessary 
one for mode I super-shear cracks in rubber-like materials [86l IIQ91 l2IQj and may well be relevant for mode II 
cracks as well. 


II for sub-Rayleigh ruptures). 


VII. FUTURE DIRECTIONS 

In the previous sections we highlighted a number of interesting topics in dynamic fracture that, in our opinion, 
have progressed significantly in recent years. Below we list a number of promising research directions that are 
related to each topic. 

Equation of motion for a crack. The work reviewed in Sect, [^has conclusively shown that the equation(s) 
of motion for simple cracks are correct, as long as a crack remains simple. The energy balance criterion 
works perfectly as long as the propagation direction of the crack is both prescribed and the relevant boundary 
conditions are properly accounted for. We do not have, as yet, a general equation for determining the path 
of dynamic cracks. One such equation was discussed in Sect. |IV C| and was applied to crack oscillations in 
2D. Other interesting work has been performed in relation to dynamic path selection in crystalline materials 
[ml EH]. Nevertheless, there is still much more to be done in this direction. Recent post-mortem analysis 
of fracture surface markings in PMMA |2I3j has suggested that microscopic crack front motion takes place 
via the merging of voids ahead of the crack tip, at velocities significantly below the mean crack front velocity. 
Whereas these observations are perfectly compatible with LEFM, they raise an interesting question of size 
effects; how large do inhomogeneities have to be in order to invalidate the continuum description of crack front 
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propagation? Whereas the interaction of crack fronts with single asperities has been explored in the LEFM 
framework [22]-[24], it will be interesting to see how/if this framework generalizes to a large number/size of 
defects. Is there a scale where the crack front loses its identity as a coherent entity? 

Nonlinear elasticity and soft solids. LEFM, as well as the weakly nonlinear theory, assume a “base” state 
of zero background strain. This approximation works well for most “standard” brittle materials, such as 
glass, brittle ceramics or brittle acrylics. Lately, new classes of materials are becoming available that join 
natural rubber as elastic, yet brittle, materials that can sustain large background strains prior to fracture. The 
fracture resistance of some of these synthetic materials, such as double-network hydrogels made of ionically and 
covalently cross-linked networks can be made to be enormous and new applications are expected to 

abound [uniiini- In addition to such “chemical gels”, there is also clear experimental evidence that elastic 
nonlinearities are important to fracture processes in “physical gels” |IQ8j . i.e. weakly cross-linked elastomers, 
that are important in biological applications. As all of these materials exhibit large strains prior to fracture, 
existing theoretical frameworks should be extended to large background strains. Very recent experiments have 
shown these effects to be important and the first extension of weakly nonlinear fracture mechanics to moderate 
background strains has been developed |222j . This framework might also apply to the related topic of the 
strength of tough pressure-sensitive adhesives, such as adhesive tapes [22311225] . 

Instabilities in dynamic fracture. The recent advances reviewed here may provide us with some tools needed 
to initiate a more complete understanding of dynamic instabilities. From these results, it appears that to 
fully understand these phenomena, the framework of fracture mechanics should be extended to systematically 
include 3D crack propagation. The near vicinity of the crack tip (e.g. process zone) is not simply a passive 
region that is needed to dissipate energy, but, as we saw in Sect. llvUl the process zone can actively feed 
back to the elastic fields driving the crack to produce unstable behavior. The 3D form of micro-branches 
implies that these effects entail a finite-size perturbation to the crack front. Initial studies of crack front 
waves considered only linear perturbations to an LEFM-described front. Further development of a nonlinear 
description of crack front behavior may provide a way to better explore crack front interactions with either 
active noise sources (e.g. micro-branches) or quenched noise, such as material inhomogeneities. Other promising 
directions to describe 3D effects include the use of phase-field models to provide an accessible way to couple 
the nonlinear elastic fields to dynamic models describing the crack tip. Some initial steps in this direction 
have been accomplished [SH 111711226] , but there are still a number of open questions about how to ensure the 
physical predictive power of these models. 

Fracture-based descriptions of friction. Recent experiments have demonstrated that frictional motion is 
intimately related to dynamic fracture. A fundamental understanding of these important processes must, 
therefore, involve rapid fracture. While many related issues has been extensively discussed in the geophysical 
literature |I42[ 1145] , careful examination of frictional “fracture” in well-controlled laboratory experiments has 
only began to emerge. In Sect. VL a few aspects of this research direction have been mentioned. Different 
modes of earthquake-like rupture states are expected to abound in these systems; these range from sustained 
slow fracture modes to super-shear rupture. The results described in Fig. [F for example, relate to dry and 
rough frictional interfaces formed by identical materials. It would be interesting to make a similar comparison 
to analytical results that describe mode II super-shear cracks along homogeneous interfaces, whose singular 
properties are predicted to be quite different from those of standard cracks [199] . Entirely different classes of 
rupture, that include “self-healing” pulse-like as well as crack-like rupture states, are believed be generated 
within bimaterial interfaces; interfaces formed by two materials having different elastic properties [11511223- 
I231j . Recent experimental results have demonstrated the existence of a number of these interesting 

modes. Crack-like modes may also be excited within lubricated interfaces [19411195] . Like super-shear states, 
these rupture modes are predicted to possess functional forms that are still singular, but with singularities that 
may be entirely different from the square root singularity generally associated with cracks. 

In conclusion, in this perspectives piece we have presented a number of new and interesting developments 
in the field of dynamic fracture. The examples noted are all important, but are by no means all-inclusive. 
We have endeavored to show that dynamic fracture is not only important in itself, but is critical to our 
fundamental understanding of topics ranging from the strength of new “super-tough” materials to the dynamics 
of earthquakes. We expect that these developments will continue into the next decades. 
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